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1.  Introduction. 

In  several  longitudinal  investigations  the  estimation  of  the  mean 
survival  time  is  of  basic  importance.  This  is  usually  based  on  the  data 
gathered  from  a  sample  of  n(2  1)  identical  units  on  test  as  in  a  clinical 
trial  or  life  test.  However,  it  is  often  the  case  in  these  survival  studies 
that  the  lifetimes  of  the  specimens  are  not  completely  observable  due  to 
random  withdrawals  and  consequent  loss  to  follow-up.  Accordingly,  for 
each  survival  time  X  we  envisage  a  competing  censoring  time  Y  ,  but  the 
only  datum  available  is  (Z,6)  where  Z  =  min(X,Y)  and  6  =  1  whenever 


X  £  Y  and  6=0  otherwise.  We  shall  assume  X.  to  be  independent  of  Y 
and  Z  &  0  almost  surely  (a.s.)  for  each  6  . 

Suppose  {(Z^,6j):l  £  i  5  n}  is  a  random  sample  of  size  n  with 
lifetimes  X^  having  common  exponential  survival  function 
F(t)  =  exp(-  t/8)  ,  t  >  0  with  0(>  0)  unknown  and  the  censoring  times 
Y^  having  common  (unknown)  survival  function  G  .  We 

A 

consider  the  sequence  of  estimators  {0^}  of  0  given  by 


i-i 


«  Ii”i  5i  2  1 

lf  h°i  -  0  • 


[l.i] 


@n  is  easily  seen  to  converge  a.s.  to  0  as  n  — for  each  6  . 


The  loss  incurred  in  estimation  of  0  by  0  is 

n 

*  2 

L  (a,c)  =  a(8  -  0)  +  cn 

n  n 


[1.2] 


where  a(>  0)  is  a  known  constant  and  c(>  0)  is  the  cost  per  unit 


observation.  With  the  definition  [1.1]  we  can  show  that 

B(0n  -  0)2  =  (aeW-1  +  o(n_1)  (1.3) 

whence  the  risk  associated  with  the  estimation  scheme  becomes 
Rn(a,c)  »  ELQ(a,c)  -  (a02/E6)n  *  +  cn  +  o(n  V  • 

Therefore,  in  order  to  minimize  this  risk  (with  respect  to  n)  we  must 
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take  a  sample  of  size  tP  with  subsequent  minimum  risk  R®  where  (for 

c  c 


small  c) 


R®  ~  2cn^  and  -*  (aB^/cES)**  . 
c  c  c 


11.4] 


Since  8  is  the  unknown  parameter  which  we  wish  to  estimate,  this  minimum 
risk  R^  and  optimal  sample  size  n^  are  also  unknown  and  thus  we  are 
led  naturally  to  explore  a  sequential  scheme  for  estimating  8  .  In  this 

article  we  propose  one  such  procedure  in  which  the  sample  size  is  given 
by  a  random  stopping  number  N(c)  and  the  associated  risk  R*  «■  EL^cj 

*  o 

satisfies  Rc/Rc  — ¥  1  as  c  — *■  0  ,  that  is  if  our  scheme  is  asymptotically 
risk-efficient. 


Sequential  procedures  analogous  to  the  one  outlined  here  have  been 
considered,  in  the  absence  of  censoring,  by  several  researchers  beginning 
with  the  pioneering  work  of  Robbins  (1959)  for  the  estimation  of  the  mean 
of  a  normal  population.  This  was  later  extended  by  Starr  (1966)  and  Starr 
and  Woodroofe  (1969) .  The  study  of  sequential  point  estimation  of  the 
exponential  mean  is  taken  up  in  Starr  and  Woodroofe  (1972) .  Gardiner  and 
Susarla  (1982)  make  the  first  examination  of  the  problem  of  sequential 
estimation  of  the  mean  survival  time  in  the  presence  of  censoring  when  both 
the  underlying  survival  time  and  censoring  distributions  are  unspecified. 
When  censoring  is  absent,  similar  procedures  for  estimation  of  functionals 
of  an  unspecified  distribution  are  discussed  in  Sen  and  Ghosh  (1981),  Ghosh 
and  Mukhopadhayay  (1978)  with  some  extensions  and  refinements  by  Chow  et. 
al.  (1981,  1982).  Throughout  the  rest  of  the  paper,  the  terms  involving 
the  random  variable  **  0]  are  left  out  without  any  further 

indication  since  ^i  “  goes  to  zero  at  an  exponential  rate  and 

all  our  scale  factors  will  be  lika  n  a  with  a  >  0  . 
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2.  The  Sequential  Procedure 

Let  m(2  1)  be  the  Initial  sample  size.  Then  motivated  by  [1.4] 
define  the  stopping  number  Nfi  (=  N (c))  by 

N  -  min{n  2  m:n  a:  (^^(o  +  n  "*))  [2.1] 

c  c  n 

«  2 

where  y  >  0  is  a  constant  and  oq  is  an  appropriate  estimator  of 
o2(0)  =  02/E6  "  03(J^  FG)  *  .  Notice  that  by  the  ordinary  Central  Limit 
Theorem  we  have  from  [1.1] 

-  8)  -*p  N(0,o2(6)).  [2.2] 

2 

An  appropriate  estimator  of  o  (8)  based  on  {CZ^,d^):l  £  i  £  n}  is 

o2  where 
n 

%  -  5X  “  >■  °)  I2-Jl 

with  the  overscore  denoting  the  usual  corresponding  sample  mean  and  [A] 

denoting  the  indicator  of  the  set  A  .  Observe  that  for  each  fixed  8  , 

«2  2 

on  converges  a.s.  to  a  (0)  as  n  , 

Our  proposed  sequential  scheme  will  take  N(c)  observations  and 

estimate  0  by  with  resulting  risk 

R*  -  a  EC8nCc)  -  6)2  +  c  E(N(c)>  .  [2.4] 

The  main  result  of  this  paper  is 

Theorem  1.  (Risk  Efficiency) .  For  each  y  >  0  , 

*  q 

R  /R  — ►  1  as  c  +  0  . 
c  c 

Furthermore,  we  can  easily  demonstrate 

a 

Theorem  2.  (Asymptotic  Normality  of  As  c  +  0  ,  for  each  8 

{ /n (c)  converges  in  distribution  to  a  standard  normal 


random  variable. 
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The  proofs  of  these  theorems  are  outlined  in  the  next  section.  He 
begin  with  a  demonstration  of  the  expansion  [1.3]  and  a  lemma  concerning 
the  rate  of  growth  of  N(c)  as  c  +  0  . 

3.  Proofs. 

He  begin  with  a  Lemma  on  left  truncated  inverse  binomial  moments. 

Lemma  1.  Let  V  be  a  binomial  random  variable  with  parameters  n  ,  p  . 
Then  for  any  positive  integers  k  ,  t 

E(V-k[V  ii])  i  K(np)"k 

where  the  constant  K  does  not  depend  on  n  and  p  . 


Proof .  Write  n^  =  n(n  -  1) . (n  -  k  +  1)  and  n-^  »  (n^)~*  .  Then 

l  {(*  +  k)'IW  .  <*>  p*(l  -  p)n~X/x! } 

x=0 

*>  (n  +  k)"Ik]  l  {(n  +  k)tx+kl  pX(l  -  p)n“X/(x  +  k)!} 
x=0 

5  Cn  +  k)”^1  p"k 

f  kl  k 

Therefore  since  (x  +  k)1  J/x  is  decreasing  with  respect  to  x  , 

E(V_k[V  a  f]  S  +  k)Ikl  l~ k  E{ Cv  +  k)"Ikl[V  *  l] } 

S  (JL  +  k)IkJ  l~k  (n  +  k)~Ik]  p'k 
S  K(np)_k  . 

We  now  turn  to  the  crucial  expansion  given  in  [1.3].  To  this  end  note  that 
from  [1.1]  we  can  write 


where 


e  -  e  -  u  v 

n  n  n 


°n  "  «n  -  e‘n>'  \  ~  K  • 


He  then  have 


[3.1] 


[3.2] 


Lemma  2.  For  any  p  k  1  ,  E J  U^J 2p  -  0(n“p)  and  E | | P  -  0(1) 
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Proof .  Since  6  is  either  0  or  1  and  0  £  Z  £  X  a.s.,  all  moments 

of  5  and  Z  are  finite.  Proa  [3.2J,  we  see  that  Dq  is  the  average 

2 

of  n  iidrv’s  with  zero  mean  and  variance  0  E5(*  0  EZ)  .  Then  by  the 
Marcinkiewicz-Zygmund-HBlder  inequality  we  obtain 
E|uJ2p  £  k^nP  e|z  -  05 1 2p  -  0(n"p)  . 

(Here,  and  throughout  this  sequel  kj  is  a  generic  constant  not  depending 
on  n) . 

For  the  second  part  of  this  lemma  we  apply  Lemma  1  to  the  binomial 

random  variable  n6Q  ,  from  which  the  required  statement  follows  immediately. 

This  concludes  the  proof  of  Lemma  2. 

Before  we  proceed  further,  we  point  out  that  [1.3]  follows  from  the 

equality  0  -  0  «*  (E5)  ^  U  {1  -  V  (6  -  E6)}  where  U  and  V  are 
n  n  n  n  n  n 

defined  by  [3.2]. 

The  next  Lemma  gives  a  convergence  rate  (as  c  +  0)  for  the  stopping 

number  N  of  [2.1].  First  note  that  N  2:  (— )1/2(1+y)  a.8.  8<J  that  n  +  » 

c  c  c  c 

a. 8.  as  c  +  0  .  In  the  sequel  int  [x]  denotes  the  greatest  integer 
£  x  and  k^.k^,***  denote  constants  independent  of  c  ,  but  could  depend 
on  G  and  0  . 

Lemma  3.  Let  0  <  e  <  1  be  arbitrary.  Then  as 
P[N  £  lnt[n°(A  -  e)]J  -  0(c) 

c  c 

and 

P[N  >  int[n®(l  +  e)]J  -  0(c)  . 
c  c 

Proof.  Write  b  -  C-)**  and  n,  -  int[b1/cl+1r)] 

c  lc 

On  the  set  [N  £  n.  ]  we  have  that  n  2  ba  for 
c  zc  n 


c  f  0 

13.3] 

»  n2c  “  lnt[n°(l  -  e)]  . 
some  n  e  (n^c. 


Therefore 


,n2c}  ‘ 
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*(NC  *  n2cJ  S  PI  U  l5n  S  bn“l11 
n*n, 
lc 

n2c 

S  P(  U  [o2  -  a2  £  b~2n?  -  a2)]  ,  where  o2  «  o2(0)  -  BZ/E6  , 

n*n,  4C 

lc 

®2c 

S  P[  0  l|c2  -  o2|  i  o2  e(2  -  e)]] 


^2  2 

«  Pl  max  |o*  -  a  \  2  nj  ,  for  some  n  >  0  .  |3. 4J 

n,  SnSn,, 
lc  2c  1 

How  a2  -  o2  “  a  (I  -  EZ)  -  b  (6  -  E6)  13.5] 

n  n  n  n  n 


where 

«n  -  (Zn  +  EZ)i~3  [7n  ^  0]  and 

bn  -  (EZ)2(E6)-3(62  +  InE«  +  (E«)2)fi  “3  lfin  +  0]  . 

-2  °  —2 

Since  a  — *■  20 (E6)  «■  a.  a.s.  and  b  —*■  SO'tES)  “  b_  a.s.  as  n  — »•  »  , 

n  ■  U  n  U 

we  obtain  from  13.5]  that  the  rlghthand  side  of  13.4]  may  be  bounded  by 

P[  max  |Z  -  EZ|  >  n]  +  P{  nax  |Tn  -  E«|  >  n]  13-6) 
nXcSnSn2c  “  nlcSnSn2c 

for  an  t|  >  0  .  Each  of  the  two  terms  in  13.6]  are  handled  In  the  same 

fashion.  Observe  that  {|Zn  -  EZ|:n^c  £  n  £  n2c^  *s  a  reverse  submartingale 

to  which  the  Kolmogorov  maximal  inequality  applies.  Thus 

PI  max  |Z  -  EZ.|  2  n]  S  k.  E|z  -  EZ|4+4y  -  0(c)  , 

n,  SnSn,  n  1  ’He 

lc  2c 

where  k^  is  a  constant  depending  only  on  y  .  Thus  combining  [3.4] 

through  [3.6]  yields  the  first  statement  of  Lemma  3.  The  second  is  proven 

in  much  the  same  way  and  so  the  proof  of  the  lemma  is  terminated. 

He  may  utilize  the  expansion  in  [3.5]  to  obtain  the  uniform  lntegrability 

of  {Nc/n^:0  <  c  £  c^}  for  some  Cq  sufficiently  small;  from  this 

E{N  /vP)  -4-1  as  c  4  0  will  follow, 
c  c 


To  this  end  note  from  [2.1}  that 


®N(c)/o  5  Nc/nc  5  °N(c)-l/o  +  (N(c)  "  1)Y/0  +  (ba)  1  * 

Then  since  N  +  •  a.s.  (as  c  +  0)  ,  we  obtain  N „/vP  — *  1  a.s., 
c  c  c 

Furthermore 

(N>°)2  s  8{(5N(c)_1/a)2  +  <f2(N(c)  -  1)"2y  +  (ba)~2)  , 

0  2 

Thus  in  order  to  show  E(  sup  (N  /n  )  )  <  «  for  some  cn  sufficiently 

-  C  C  v 

CSCq 

^2  2 

small,  we  need  only  examine  E(sup|S  /a  -  1)  for  sufficiently  large  m  . 

n>m  n 

However,  from  [3.5]  and  the  convergences  afl  — ►  a^  ,  bQ  — *■  b^  a.s.  it 

suffices  to  establish  the  finiteness  of  E(sup|Z  -  Ez|)  and 

n£m 

E(sup|6  -  EdJ)  .  The  second  is  trivial,  while  the  first  is  bounded  by  four  times 
nStra  ” 

_  2  _ 

Elz  -  Ezl  since  {lz  -  EZl  :n  S  m}  is  a  reverse  submartingale. 

*  •  *  n  * 

Proof  of  Theorem  2.  Recall  [3.1]  and  [3.2].  From  the  central  limit  theorem 


n*Un  — *-p  N(0,82E5)  ,  for  each  8  >  0  . 


Also  V  — *■  (E6)  a.s.  Hence 
n 


i  ^  2  *  ? 

n*(8  -  8)  — N(0,o  )  .  We  have  noted  that  a  — *■  a  a.s.  (as  n  — *■  <*•) 
n  V  n 

and  N£  f  as  c  1  0  .  Hence,  by1  [3,8],  [3,13],  and  [3,15J  belpw,  we  get 


(0N(c)  -  0)/oN(c)  -+V  »<0,1>  *s  c  +  0 


Proof  of  Theorem  1.  In  view  of  the  fact  that  E{N_/n  }  — *■  1  ,  [1.4]  and 
-  c  c 


[2.4]  we  must  verify  that 

lim(cn®)-1  E(§n^cj  -  8) 2  ■  1  . 

Recall  the  expansion  [3.3].  We  first  show 

lim  (cn®)’1  E((8n(c)  -  8)2[Nc  S  n2c]}  -  0 


[3.7] 


[3.8] 


E{(eN(c)  -  6)2[Nc  S  n2c]>  *  k3  [EU2(c)[Nc  S  n2c]  +  EW2(c)[Nc  £  n2c])  ,  [3.9| 
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Now  applying  the  maximal  inequality  to  the  reverse  submartingale 
{IV  eTJ:nlc  *  n  5  a2c}  in  the  third  step  and  Lemma  3  we  have 

E{UH(c)|Nc  S  V”  S  E(  *“  »>cS”2c» 

"lc^n2c 


£  gl/(l+Y)(  max  |jj  |2+2Y)(P[N  ' -i  n, J)y/(1+Y> 
nlc^a2c 

5k.E1^l+Y^|u2+2Y|  0  (cY/(1+Y>) 


‘lc 


0(n~J)0(cy/1+Y)  . 


Since  nlc  ~  c-1/2(1+y)  and  (cn®) ~  c-ls 


we  have  that 


Um(cu»rl  E(uJ(c)[Hc  S  n2c)}  -  0  . 

n2c 

Similarly  E{W2(c)[Nc  £  n^]}  -  l  E{W2[Nc  -  n]} 

nanlc 


[3.10] 


2c 


s  (  I  E  W2P)1/PCPIN  5  n«  ]) 


1-1/p 


n=n 


2c J 


P  >  1 


[3.11] 


lc 


Now  E  W21*  -  E{U2p  V2p(6  -  E6)2p} 

n  n  n  '  n 

S  E1/3{U6p}E1/3{V6p}E1/3(6  -  E5)6p 

n  n  n 

-  {0(n~3p)0(l)0(n~3p)}1/3  =  0(n~2p)  , 

Thus  the  rhs  in  [3.11]  is  0(n"2+1/p  c1-1/p)  =  0(ch)  ,  where 


h  »  {2(1  +  y)}  *(2  -  1/p)  +  (1  -  1/p)  >  whenever  p  >  (2y  +  3)/(3  +  y)  .  So 

[3.12] 

Collecting  our  results  [3.10]  and  [3.12]  we  obtain  [3.8].  An  entirely 


UmCcn")-1  E«J(c)  [Hc  S  „2c)}  -  0 


analogous  argument  using  the  second  part  of  Lemma  3  will  show 

f2, 

c+0 


lim(cn^)-1  E{eN(c)  -  8)*[Nc  >  n3(J>  -  0  ,  [3.13] 


where  n3c  ■  lnc(l  +  e)]  .  Therefore,  in  order  to  verify  [3.7]  we  are 
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left  with  proving 


IWo®)'1  EU5n(c)  -  »)2[n2c  <  Nc  S  »3c]>  -  0  . 


13.14) 


0  -i  *  2 

However,  from  [1.3)  lim(cn  )  E(0  -  -  0)  ■  1  ,  and  further  from  the 

c+0  c  n° 

£ 

inequality 


El<Vc)  -  6  0><6  0  -  e)I”2c  <  »c  4  ”3c)} 
n  n 
c  c 


S  E%{<Vc>  '  lo)2[r,2c  *  »c  5  “3c»  0  -  e>2  • 


we  see  that  it  suffices  to  establish 


W<5h(o,  -  •  0)2ln2c  <  Nc  S  n3c])  -  0 

c 


13. 


Once  again  we  utilize  the  expression  [3.3]  from  which  we  get 


.-1, 


Vc)  -  9„0  -  <Vc>  -  V  -  %(c)  -  «  0» 

c  •  nc  c 


[3.16] 


Observe  that 


E«eH(c>  ‘  \0>  [n2c  *  “c  S  "3c11 
c 


S  k5'E{(U„(c)  -  0  0>  [n2c  <  Sc  5  »3c» 
v  n 

c 


[3.17] 


Now  ' 


+  MVc)  -  *  0>  t"2c  <  »c  5  "3c,)J 
c 


E«UN(c>  -  U  0>  ln2c  <  .Nc  *  n3c» 
nc 


[3.18] 


5  max  [E{  max  (U  -  U  Q)2}  ,  E{  _  max  (U  -  U  Q)2}) 

n0  <n£n  n  n  <n:£n,  nu 

2c  c  c  c  3c  c 


Since  {U:n-  <  n  £  n  }  is  a  reverse  martingale,  {(U  -  U  rt)  :n„  <  n  S  n  } 

u  4C  C  fi  U  4C  C 


is  a  reverse  submartingale  to  which  the  maximal  inequality  will  be  applied. 


This  yields 


E{ 

n 


“*  0(0n  -  ”  0>2>  5  4  E1V 

_  <nsn  n  2c 

2c  c  c 


*  4  {EU  2  -  EU  2}  13.19] 

n2c  n" 

2  2  *“1  o  a2 

Now  EUn  =  (0  E6)n  ,  n2c  =  int[n°(l  G>]  and  cn°  =  0(1)  .  Therefore  the 
rhs  of  [3.19]  is  of  order  (cn^eCl  -  e)  ^OCl)  .  Hence,  by  choosing  e 
arbitrarily  small  we  have  that 

lim(cn®)_1E{  max  (U  -  U  n)2}  =  0  .  [3.20] 

c+0  n_  <n£nu  nu 

2c  c  c 


The  second  term  in  [3.18]  is  handled  in  an  entirely  analogous  fashion 
and  will  yield  the  result  paralleling  [3.20].  So 


Finally 

handled 


limCcn")-1  E«Ub(c)  -  iy2tV  <  »c  s  «3cJ)  -  0  .  [3.21] 

c 

we  must  dispose  of  the  second  term  in  [3.17].  This  is  easily 
2  -2 

using  the  bound  EW^  =  0(n  )  and  the  simple  inequalities 

E((U»(c)  -  ”  <A2c  <  Ne  S  "3c!) 

n 

c 


n, 


3c  -2  -1 

s  4  {EW2  +  T  EW2}  =  0(n°  )  +  0(n°  )e(1  -  e)-2  . 
0  “  n  c  c 


n  n=n„ 
c  2c 


Then  on  selecting  c  arbitrarily  small  we  obtain 

limCcnJ)'1  E{(Wb(c)  -  wy2[nZc  <  Nc  S  „3c))  -  0  .  (3.22] 

c 

Therefore  in  conjunction  with  [3.17],  [3.21],  and  [3.22]  yield  [3.15]. 
This  concludes  the  proof  of  Theorem  1. 
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4.  Concluding  Remarks. 

To  start  with,  it  is  worth  pointing  out  that  the  results  are  true 
even  if  G  is  degenerate  at  any  point  T  .  That  is,  the  results  hold 
even  in  the  fixed  point  truncation  of  .  Although  the  results  are 
stated  here  only  when  are  i-i-d.,  similar  results  hold 

even  if  Y^,***,Yn,»»*  are  assumed  to  be  independent  only.  In  the 


#00  — > 

latter  case  however,  some  sort  of  condition  like  lim  Jq  FGr  >  0,  with 


n  G  «*  J  G  ,  needs  to  be  imposed. 
n  i=l  1 


n-x» 


— y 

The  final  comment  concerns  the  damping  factor  n  introduced  in  the 


definition  (2.1)  of  the  stopping  random  variable  N£  .  Without  this  factor. 


the  sequential  estimator  0^  may  not  be  asymptotically  risk  efficient 

c 


as  the  following  argument  suggests.  Define 

N 


•»c  i  55  min{n  i  m:n  2:  (-)**  on> 


[4.1) 


where  On  is  defined  by  [2.3).  With  these  definitions,  it  can  be  seen 


that  if  the  censoring  distribution  G  has  positive  mass  pat  zero,  then 


*  *  2  m 

the  risk  R£  ^  of  the  sequential  procedure  is  at  least  0  p 


c,l 


and  R  *  /R^  — *•  »  as  c  +  0  and  hence  0  is  asymptotically  totally 

Ca«L  C  W  1 

C,1 


risk  inefficient.  Thus  the  factor  n  ^  in  the  definition  of  N£  appears 


to  be  important. 
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